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1.0  INTRODUCTION 

The  Dynamic  Economic  Values  Model  (DYNEVAL)  uses  Lagrange  Dynamic 
Programming  to  chart  through  time  the  optimal  trajectory  of  an  economy 
which  has  either  been  disrupted  in  some  way,  or  in  which  people  have 
changed  their  value  structure  (e.g.,  by  giving  more  emphasis  to 
leisure  or  the  military).  "Optimal  trajectory"  is  that  time-dependent 
mix  of  investment  and  production  activities  which  allows  consumers 
(people,  government,  military,  etc.)  to  maximize  over  time  a  sum  of 
value  functions.  Too  much  early  investment  causes  people  in  early  time 
periods  to  have  too  low  a  standard  of  living  based  on  their  value 
functions;  too  little  investment  causes  the  economy  to  recover  too 
slowly. 

The  basic  logic  flow  of  the  model  is  shown  in  Fig.  1-1.  Since 
each  problem  the  model  solves  is  different,  a  storage  allocator  is 
used  to  parse  a  large  storage  array  into  a  series  of  smaller  data  arrays 
of  appropriate  dimensions  for  the  problem.  The  basic  data  for  the  prob¬ 
lem,  including  current  activity  levels  and  a  direct  requirements  matrix, 
are  then  input  to  the  program  which  proceeds  to  calculate  both  the 
current  day  economy  and  the  final  equilibrium  economy.  In  some  studies, 
one  does  not  care  about  the  trajectory  to  equilibrium,  but  rather  only 
about  the  final  equilibrium  economy  itself.  In  these  cases,  the  model 
may  be  stopped  here. 

Once  the  above  is  accomplished,  the  model  starts  its  real  job  of 
calculating  the  optimal  investment,  production,  and  consumption  activity 
levels  which  gradually  move  the  economy  to  its  equilibrium  state.  This 
is  done  by  optimizing  each  time  period  in  turn  and  propagating  the  capital 
resources  produced  forward  in  time,  taking  into  account  population  growth, 
gestation  time,  depreciation,  etc.  Each  time  period  is  presented  with 
a  set  of  capital  which  the  past  has  made  available  to  it,  and  a  set  of 
capital  values  which  the  future  is  willing  to  pay.  From  these,  and 
internal  requirements  detailed  later,  the  single  time  period  produces. 
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Figure  1-1.  Logic  Flow  of  Economic  Model 
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invests,  and  consumes  in  an  optimal  way.  As  each  time  period  considers 
its  options  it  automatically  assigns  a  rental  (shadow)  value  to  the 
capital  and  labor  needed  to  produce  the  goods  and  services  for  this 
time  period.  These  rental  values  can  then  be  propagated  backward  in 
time  (taking  into  account  discount  values,  population  growth,  gestation 
times,  etc.)  to  give  each  time  period  a  more  realistic  view  of  the 
future.  Each  time  period  is  then  reoptimized  taking  into  account  this 
new  information  and  new  capital  passed  to  it  by  the  past  and  the 
process  is  repeated. 

The  following  sections  give  the  mathematical  details  of  each  stage 
of  the  program.  These  are  technical  in  nature  and  intended  for  use  by 
someone  interested  in  understanding  the  precise  way  the  model  operates 
or  interested  in  making  changes  to  the  model.  For  the  person  interested 
in  using  the  model  as  a  research  tool.  Appendix  A  is  a  user's  manual 
explaining  how  to  run  the  model  and  giving  examples  of  output  appearing 
on  the  terminal  screen  and  the  line  printer.  Appendix  B  is  a  description 
of  the  companion  program,  AGGRAT,  used  to  aggregate  economic  data  and  put 
it  into  model  format.  Definitions  of  all  pertinent  variables  and  arrays 
are  given  in  Appendix  C.  Appendix  D  is  a  complete  listing  of  the  source 
code  for  the  model . 
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2.0  FUNDAMENTAL  RELATIONS 

Throughout  the  course  of  the  program  certain  variables  are  known 
and  others  must  be  calculated.  In  the  single  time  period  optimizer,  e.g. 
the  amount  of  capital  resources  available  for  production  is  known.  In 
the  equilibrium  calculation  the  production  and  investment  activities  are 
known  as  a  simple  function  of  the  resources,  and  the  rental  value  for 
capital  is  known  as  a  simple  function  of  its  cost.  In  order  to  calculate 
missing  quantities  certain  fundamental  relationships  among  resource 
levels,  activity  levels,  costs,  prices,  and  values  need  to  be  derived. 
This  section  details  relationships  among  these  variables  for  the  pro¬ 
duction  and  consumption  activities.  In  addition,  other  relations  are 
valid  when  the  economy  is  in  equilibrium.  These  will  be  discussed  in 
Sec.  5.  Investment  and  trade  activities  present  special  problems 
and  will  be  discussed  in  Sec.  6. 

2.1  FUNDAMENTAL  VARIABLES 

The  following  is  a  list  of  the  fundamental  variables.  Other 
subsidiary  variables  will  be  defined  as  they  occur.  Superscripts  refer 
to  types  of  activities,  i.e.,  production,  investment  (production  of 
capital),  or  consumption,  and  are  not  dummy  indices. 

K,  *  Inventory  level  of  capital  resource  j.  These 

J 

are  physical  items  (machines,  etc.)  and  may  be 
used  only  in  the  production  of  a  single  consum¬ 
able 

pff 

=  "Effective  level"  (explained  later)  of  the 

total  capital  mix  used  to  produce  consumable 

i.  To  reduce  cumbersome  notation  the  "eff" 

will  sometimes  be  omitted  in  this  text.  The 

eff 

subscript  i  will  implicitly  refer  to  ,  while 
the  subscript  j  will  refer  to  a  particular  capital 
resource  «• 

J 

L. j  =  Amount  of  labor  used  in  the  production  of 

consumable  i 


X? 


*{ 


Total  labor  pool 

Activity  level  for  production  of  consumable  i 

Activity  level  for  the  ith  consumption  activity. 

This  activity  need  not  consume  any  of  consumable 
i,  nor  is  it  restricted  to  consuming  only  one 
consumable.  Rather  it  is  defined  as  the  consumption 
of  any  particular  group  of  consumables  by  any  or 
all  of  a  group  of  consumers  (people,  government, 
etc. ). 

Activity  level  for  investment,  i.e.,  the  production 
of  capital  resource  j 

Value  of  the  effective  capital  needed  to  produce 
consumable  i 

Value  of  capital  resource  j 

Price  of  consumable  i 
Rental  value  of  capital  resource  K. 

J 

=  Rental  value  of  labor.  It  has  the  same  value 
regardless  of  the  activity 

There  are  three  basic  types  of  activities:  production  (X?), 
consumption  (X^,  and  investment  (Xt).  The  production  activities 
require  capital,  labor,  and  a  specified  mix  of  various  consumables. 

Each  production  activity  produces  only  one  tyKe  of  consumable  and  each 
consumable  is  produced  by  only  one  type  of  production  activity.  Con¬ 
sumption  activities  require  no  capital  and  no  labor  (except  consumption 
of  leisure  which  requires  only  labor),  but  do  require  a  specified  mix 
of  consumables.  The  same  consumable  may  be  used  in  several  different 
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consumption  activities.  I  v_ounent  activities  produce  capital.  They 
require  no  labor  but  do  require  a  specified  mix  of  consumables.  The 
investment  activities  are,  in  effect,  bookkeeping  activities  whereby 
various  consumables  are  turned  into  a  single  capital  resource  specific 
to  a  given  production  activity.  The  labor  required  to  produce  the  capital 
is  implicit  in  the  labor  required  to  produce  the  consumables  wed  to 
create  the  capital. 

2.2  BASIC  ASSUMPTIONS 

In  order  to  calculate  any  relations  among  the  above  listed  variables, 
certain  assumptions  must  be  made  about  how  the  economy  responds  to  changes 
in  prices,  supply,  or  demand.  The  following  assumptions  are  basic  to 
the  model.  In  some  cases  these  assumptions  are  given  functional  form. 

Where  the  form  and  its  implications  are  spread  throughout  the  model, 
as  opposed  to  a  localized  function  which  could  easily  be  changed  in  , 
the  model,  it  will  be  noted  as  such. 

1.  Consumption  activity  levels  vary  inversely  with  cost. 

The  higher  the  cost  of  the  consuming  activity,  the  lower  the  level 
of  the  activity.  If  consumption  activity  i  uses  ai •  of  consumaole 
j  per  unit  activity,  then  its  cost,  C. ,  is  given  as: 


ci  *  Z  °u  H 

j 

In  the  model  Assumption  1  is  given  functional  form  by  defining  the 
consumption  activity  level  as: 


where  F..  =  normalization  constant  (based  on  current  economy)  and  b^  = 
elasticity  of  demand  factor. 
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The  exact  form  of  Eq.  2-2  is  not  required  by  the  model.  However, 
this  form  is  a  particularly  simple  yet  useful  way  of  expressing  Assump¬ 
tion  1.  It  is  equivalent  to  defining  a  utility  for  consumption  activity 
i  of  the  form 


u(xj) 


ri 


c(1'bi) 

xP  1  -  1 

rTTT 


ri  In  (x‘)  , 


.  bi  M 

(2-3a) 

bi  =  1 

(2-3b) 

and  stating  that  xP  is  the  level  where  the  marginal  utility,  dU/dx9, 
is  equal  to  the  cost.  Usually  is  set  equal  to  1.0.  For  this  value 
of  b.j ,  doubling  the  activity  level  results  in  the  same  increment  of 
utility  regardless  of  the  original  value.  It  states,  for  example,  that 
two  people  whose  salaries  are  doubled  experience  the  same  increase  in 
happiness  regardless  of  their  original  salaries.  Values  of  b^  between 
0.5  and  2.0  seem  to  be  most  reasonable. 


2.  The  production  of  consumable  i  may  be  considered  a  function 
of  capital  and  labor  alone:  xP  =  f(kPT  ,  ) .  Furthermore,  capital 

and  labor  may  be  traded  for  one  another  in  the  production  of  any  consum¬ 
able  according  to  a  particular  relation  known  as  the  translog  function 
(also  called  the  "constant  elasticity  function"): 


v  K;si 

Ki  1 


(2-4) 


where  g^  =  normalization  constant  (based  on  current 

economy) 

yK.  +  yL.  =  1.0 

o — j-  =  elasticity  of  substitution 

»i  1 

The  implications  of  the  translog  function,  e.g.,  its  homogeneity 
[xP(aK,  aL)  =  axP(K,L)],  are  spread  throughout  the  model.  It  would  be 
difficult  to  change  to  another  production  function. 
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3.  Capital  is  used  in  place  of  labor  where  capital  is  cheaper, 
and  labor  is  used  in  place  of  capital  where  labor  is  cheaper. 
That  is,  capital  and  labor  are  traded  one  for  the  other  until 
the  value  of  their  marginal  products  are  equal.  This 
assumption  is  expressed  mathematically  as: 


3X?/3K®ff 
L  3X?/3L. 


*  constant 


(2-5) 


4.  The  price  of  any  consumable  is  equal  to  the  average  unit 
cost  of  producing  it,  i.e., 

(2-6) 


where  aij  =  flxeC*  amount  consumable  j  required  to  produce 

one  unit  of  consumable  i. 

While  capital  and  labor  are  traded  one  for  the  other  in  production 
activities,  there  exists  no  mechanism  in  the  model  for  substituting 
one  consumable  for  another,  i.e.,  the-a-j  in  Eq.  2-6  are  fixed  for  the 
given  problem. 

2.3  DERIVED  RELATIONS 

The  following  relationships  are  derived  from  the  basic  assumptions. 
Their  use  in  the  model  will  be  discussed  in  later  sections. 

1.  Assumption  1  allows  an  independent  calculation  of  A|_,  the 

rental  value  of  labor.  For  the  form  of  Assumption  1  described 
by  Eq.  2-2, 
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i 

i 

I 

i 


r 

X  =  Leisure 

L  (L  -  I  L.  ^Leisure 
i 

since  the  cost  of  leisure  is  equal  to  the  cost  of  labor. 


2.  Differentiating  the  translog  function  of  Assumption  2 

(Eq.  2-4)  and  substituting  it  in  Eq.  2-5  of  Assumption  3 
yields  a  relation  between  the  ratio  of  capital  to  labor 
and  the  ratio  of  their  rental  values: 


yK.XL 

Vs 


(2-8) 


3.  The  fair  price  assumption.  Assumption  4,  may  be  written 
as: 


Vk, 


tLi\ 


where 


a;  =  Ai 


aijAj 


(2-9) 


stating  that  value  added  is  divided  between  capital  and 
labor.  Equations  2-4  and  2-8  may  be  substituted  in 
Eq.  2-9  to  yield: 


(2- 10a) 
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3.0  REQUIRED  INPUT  DATA 

The  model  expects  the  following  data  to  be  input  at  the  start  of 
the  model.  The  first  five  items  describe  the  given  economy  directly 
and  cannot  be  changed  within  the  model.  The  last  five  are  particular  to 
the  study  at  hand  and  may  be  changed  for  each  run. 

1.  Current  day  activity  levels.  These  include  the  investment, 
trade,  and  consumption  activities  as  well  as  the  production 
activities. 

2.  Current  day  capital  resource  inventories. 

3.  Gestation  times,  depreciation  rates,  discount  rates,  and 
growth  rates  (above  population  growth)  for  all  capital 
resources.  These  values  are  usually  calculated  or  input 
at  the  time  the  model  data  is  formed  by  using  the  date 
aggregator  AGGRAT,  which  is  described  in  Appendix  B. 

4.  A  direct  requirements  matrix,  a.  giving  the  amount  of 

•  w 

consumable  j  used  per  unit  of  activity  i  for  all  activities 
i  mentioned  in  (1)  above. 

5.  A  rental  array  giving  the  current  mix  of  capital  and  labor 
for  each  of  the  production  activities. 

6.  The  values  of  b^  used  in  the  utility  function  for  consump¬ 
tion  activities  (see  Assumption  1*.  Sec.  2). 

7.  The  fraction  of  each  capital  resource  which  "survived" 
the  proposed  disruption. 

8.  Minimum  consumption  activity  levels.  While  this  array  is 
normally  set  to  zero,  it  allows  the  user  to  specify,  e.g., 
that  a  minimum  of  agricultural  products  must  be  consumed. 

By  setting  the  minimum  activity  level  negative,  the  user 
allows  a  zero  level  of  consumption. 

9.  The  ratios  of  two  utility  coefficients  T;  those  used  for  the 
study  to  the  current  day  r  values.  If  no  change  in  priorities 
is  postulated,  these  ratios  are  unity. 
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The  translog  elasticity  (see  Sec.  4.2)  parameters 
8  for  both  the  capital  labor  trade-offs  and,  if  two 
types  of  capital  are  required,  the  effective  capital 
calculation. 
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4.0  CURRENT  ECONOMY 

While  the  inventory  and  activity  levels  of  the  current  economy  are 
input  directly  to  the  model,  the  utility  function  and  translog  coeffi¬ 
cients  which  describe  the  economy  must  be  calculated.  This  is  done  by 
subroutine  CURRNT.  The  logical  flow  of  CURRNT  is  diagrammed  in 
Fig  .  4-1. 


4.1  UTILITY  COEFFICIENTS 

Equation  2-2  gave  the  consumption  activity  level  as  a  function  of 
the  cost  for  that  activity. 


The  Ti  are  normalized  to  the  current  economy.  Since  physical 
units  are  defined  as  the  amount  valued  at  1.0  in  the  current  economy, 
Ci  *  1,  and 


r 


i 


(4-2) 


where 


and 


Z.  =  minimum  activity  level  (offset), 

F^  =  ratio  of  study  economy  to  current  economy 
priorities, 

Q 

Xi  refers  here  to  the  current  consumption  level. 


The  parameters  b^  are  input  by  the  user  (see  Appendix  A). 


4.2  TRANSLOG  PARAMETERS 

Capital  and  labor  may  be  traded  one  for  the  other  according  to 
the  translog  function,  Eq.  4-4: 


9i 


0 


n 


(4-3) 


Figure  4-1.  Logic  Flow  of  Current  Economy  Parameter  Calculate 


The  elasticity  parameter,  8..,  is  input  directly  or  indirectly  by 
the  user  (see  4.2.3  below).  The  y's  and  e^  must  be  calculated. 

4.2.1  Capital  and  Labor  Coefficients  y 

Equation  2-8  gave  the  relation  between  the  ratio  of  capital  and 
labor  and  the  ratio  of  their  rental  values  as: 

S.+l 


ft) 


=  Yk1xL 

\\ 


(4-4) 


In  the  current  economy,  X^  =  1.0  by  definition  of  the  units  of 
labor.  The  rental  value  for  capital  may  be  found  by  dividing  the  total 
spent  by  industry  i  on  capital  by  the  total  inventory  of  the  capital 
used  by  the  industry: 


(4-5) 


where 


r  *  amount  spent  on  capital  j  in  the 
current  economy. 


Since  K.  and  L.  are  known  for  the  current  economy,  as  well  as 

II  l\* 

and  X,  ,  they,,  and  y.  may  be  found  from  Eq.  4-4  and  the  fact  that 
Ki  Li 

their  sum  must  be  unity: 


6,+l 

x*  (K./U  ’ 


Ki '  i  i 


si+l 

xK^Ki/Li^  +  \ 


(4-6a) 


(4-6b) 


15 


4.2.2  Normalization  Factor  g.. 

The  parameter  g^  is  simply  the  value  required  to  reproduce  the 
current  level  of  activity  if  the  capital  and  labor  used  equal  that  used 
in  the  current  economy. 


9  . 


l 


X? 


(4-7) 


4.2.3  Elasticity  Parameter  ei 

Under  most  circumstances  the  elasticity  parameter  is  input  by  the 
user.  Valid  values  for  8^  lie  in  the  range  -1.0  £  8^  £  ».  If  8,.  lies 
outside  this  range  the  model  assumes  the  user  wants  the  model  to  calcu¬ 
late  the  actual  6^  and  that  the  input  value  refers  to  the  productivity 
at  infinite  labor  levels.  Some  industries  are  able  to  increase  their 
output  significantly  by  increasing  the  amount  of  labor  used,  e.g.,  by 
adding  a  second  shift.  Other  industries  are  not.  If  the  user  thinks 
that  the  most  a  particular  industry  could  increase  its  output  given  fixed 
capital  and  unlimited  labor  is  by  a  factor  p  >  1.0,  he  may  set  8^  =  -p, 
and  the  model  will  calculate  8^  such  that: 


(4-8) 


Equation  4-8  may  be  rewritten  as: 


(4-9) 
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From  Eq.  4-4: 


(4-10) 


or 

1n  yyy 

"ln  ak(Ki/Li ^  +  1,0 
Si  In  p 


(4-11) 


4.2.4  Effective  Capital 

Normally  a  given  production  activity  uses  one  type  of  capital. 

The  model  allows  two  types  of  capital  per  industry,  however,  and  treats 
the  combination  as  an  "effective"  level  of  capital.  The  two  types  of 
capital  are  combined  by  a  separate  translog  function: 


(4-12) 


To  avoid  cumbersome  notation,  we  will  set  and  "j^"  and  "J2" 

to  "1"  and  "2"  respectively  in  this  section. 
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Eq.  2-5  may  be  generalized  for  this  case  to  derive  the  relationship 
between  AK^and  the  individual  rental  values: 


where  j  =  1  or  2. 


(4-13) 


While  Eq.  4-13  may  look  imposing,  it  offers  a  very  simple  repre 

sentation  for  A„  given  A,,  and  A„  : 

*eff  N  *2 


XKxK1  +  XK2K2  =  e(Keff/e^  1  (YiKl  2  +  y2  ^XK. 


■eff 


KeffXKeff' 


or. 


X1K1  +  X2K2 
Keff 


(4-14) 
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Equation  4-14  states  that  the  effective  rental  value  of  the  capital 
combination  is  equal  to  the  total  spent  on  both  divided  by  the  effective 
level  of  capital.  This  property  follows  directly  from  the  homogeneity 
of  the  translog  production  function. 

The  model  calculates  the  y  and  g  parameters  for  the  effective 
capital  function  in  a  way  similar  to  Secs.  4.2.1  and  4.2.2.  Equation  4-5 
gives  the  ^Kj  for  the  current  economy.  By  replacing  with  K^,  with 
K2’  ^Ki  ^i  and  \  ^2’  a^ow  the  calculation  of  y^ 

and  Yk  ..  Similarly,  replacing  X?  with  ^  and  B.j  with  B,  Eq.  4-7 
allows  the  calculation  of  g. 
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5.0  EQUILIBRIUM  ECONOMY 

The  equilibrium  economy  is  defined  as  that  economy  where  capital 
is  being  produced  at  such  a  rate  that  capital  per  unit  population  re¬ 
mains  constant.  Furthermore,  real  prices  must  remain  constant.  For 
these  reasons,  the  equilibrium  economy  is  not  necessarily  the  current 
economy,  although  it  could  be  if  no  industry  is  experiencing  a  growth 
greater  than  population  growth  and  there  are  no  consumer  preference 
changes.  The  equilibrium  economy  is  calculated  in  subroutine  EQLBRM. 
The  logical  flow  of  EQLBRM  is  diagrammed  in  Fig.  5-1. 


5.1  CAPITAL  INVENTORY  LEVELS  AND  INVESTMENT 

Capital  must  increase  with  the  population  to  provide  the  same 
level  of  production  per  person.  Assuming  an  exponential  growth  in 
population,  P,  and  a  depreciation  rate,  d,  the  capital  inventory  at 
time  t,  K(t),  is  related  to  the  inventory  at  time  t  by: 


K(t)  =  K(t0 


)  e 


P(t 


■v 


(5-la) 


and 


-d(t-t  ) 

K(t)  =  K(t0)  e  0  + 


-  t)  e 


-d(t-t') 


dt ' 

(5-lb) 


where  t  =  gestation  time 

X  (t1  -t)  =  amount  of  new  capital  becoming  available  at  t' 

Note  that  Eq.  5-lb  is  a  statement  of  physical  fact.  Inventory 
at  time  t  must  be  equal  to  the  inventory  at  time  t  which  has  not  depre¬ 
ciated,  plus  what  has  been  made  between  t  and  t.  Eq.  5-la,  on  the 
other  hand,  is  valid  only  in  equilibrium.  It  states  that  capital 
inventory  increases  at  the  same  rate  as  population. 
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Figure  5-1.  Logic  Flow  of  Equilibrium  Calculation 
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Setting  the  right  hand  sides  of  Eqs.  5-1  equal  to  each  other  and 
differentiating  with  respect  to  t  yields  the  expression: 


K( 


f  P(t-t  )  -dU-Ol 

t„>  [pe  +de  0  J- 


X!( t  -  x)  -  d  /V(f  -  t)  e'^*'*' 


dt ' 


(5-2) 


The  integral  in  Eq.  5-2  may  be  related  to  K(t)  and  K(tQ)  from  Eq.  5-1 
to  finally  yield: 


K(tQ)(P  +  d)  e 


„J 


X‘(t  -  t) 


or 


P(t-tJ  Pt  T 
K(tJ(P  +  d)  e  0  ePT  =  X1 


(5-3) 


(t) 


which,  through  Eq.  5-la  becomes 

K(t)(P  +  d)  ePT  =  X1 (t)  (5-4) 

Equation  5-4  relates  the  activity  level  in  equilibrium  to  the 
inventory  level.  The  inventory  level  will  grow  in  time  at  the  same 
rate  as  population.  It  is  convenient  in  the  model  and  clearer  to  the 
user  to  normalize  inventory  and  activity  level  to  the  current  day 
population  so  that  the  time  dependence  in  Eq.  5-4  may  be  removed. 

After  all,  there  is  no  way  to  know  before  the  model  is  run  how  long  it 
will  take  to  reach  equilibrium,  and  a  curve  (trajectory)  which  flattens 
at  equilibrium  is  much  easier  to  interpret  than  one  which  grows 
exponentially  at  the  population  growth  rate.  With  this  normalization, 
Eqs.  5-1  become: 
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K(t)  =  K(t0)  (5-5a) 

K(t)  =  K(t  )  e'(P+d)(t'to)  ♦  e-PT  f%-  -  t )  e-<P+D)ft-f)  dt. 

-'to 

(5-5b) 

where  K  is  normalized  to  population  and  is  not  the  same  as  the  K  of 
Eqs.  5-1.  The  end  result  of  solving  Eq.  5-5  for  X*  as  a  function  of  K 
is  identical  to  Eq.  5-4  with  the  time  dependence  removed: 

*  Equilibrium  “  ^Equilibrium^  +  (5-6) 

5.2  CAPITAL  RENTAL  VALUES 

In  equilibrium,  the  rental  value  for  capital  is  related  to  its 
cost  of  production  in  a  very  simple  way.  Since  the  cost  of  production 
must  be  equal  to  the  sum  over  time  of  all  the  discounted  rents,  and 

since  in  equilibrium  the  rental  value  must  remain  constant,  Eq.  7-10 

(See  Sec.  7)  reduces  to: 

A  =  X  e-(P+p)T 

Equilibrium  Equilibrium  p--+  ^  +-—  (5-7) 

where  A  =  cost  of  production 

X  =  rental  value  per  unit  time 
p  =  discount  rate 

5.3  FINDING  THE  EQUILIBRIUM  STATE 

The  following  is  the  algorithm  used  by. the  model  to  calculate 
the  equilibrium  state  (see  Fig.  5-1).  It  is  an  iterative  technique 
which  converges  rapidly.  If  the  curr w.  economy  contains  no  growth  relative 
to  population  growth,  i.e.,  is  itself  an  equilibrium,  and  if  the  utility 
parameters  have  not  changed,  the  process  clearly  takes  one  iteration. 

If  the  above  is  not  true,  the  process  may  take  ten  to  twenty  iterations. 


23 


The  reader  should  note  that  by  this  time  the  parameters  y,  3,  r, 
b,  and  a  have  been  calculated  from  the  current  economy  and  are  fixed 
hereon  and  available  to  the  model. 


1.  Assume  some  set  A^  of  consumable  values  and  some  set  K.  of 
capital  resources,  e.g.,  the  current  economy  values. 

2.  Calculate  the  investment  activities  xj  from  Eq.  5-6. 

3.  Find  the  cost  of  capital,  a|,  from  the  consumable  values  A. 

4.  Calculate  the  rental  value  X,,  ,  from  A?  and  Eq.  5-7. 

Ki  1 

5.  Calculate  the  labor  required  for  each  production  activity 
from  Eq.  2- 10a: 


(5-8) 


where  Al  =  A.  -  I  a..  A. 

j  1 J  J 


and  a.  .  =  amount  of  consumable  j  required  per  unit  activity  i 
1  j 


6.  From  L.  and  K..  calculate  the  production  activities  x!j  from 
the  translog  function,  Eq.  2-4: 


P  ~pi  "pi 

Xi=ei  *K*i 


(5-9) 


7.  Find  X^  from  the  utility  function  for  leisure: 


Leisure _ 

~  ^Leisure 


(L  -E  L.) 


(5-10) 


8.  Find  a  new  estimate  for  from  Eq.  2-6: 


*  Calculated  _ 
i 


X  K  +  XL  +  z  X^a-.A. 
=  i  i  jyH  1  J 

0  -  af1)  X? 


(5-11) 


New  .  Calculated  .  . 
A"eW  ■  H  (A •  +  Ai ) 


(5-llb) 


9.  If  Aj  =  A^,  proceed.  Otherwise  go  back  to  Step  3. 

10.  Next  find  the  cost  per  unit  consumption  activity,  and  the 
level  of  that  activity: 


Ai  ■  j  «ij*j  . 


(5-12a) 


(5-1 2b) 


11. 


p 

Compare  how  much  of  each  commodjty  was  produced,  x^ ,  with  how 
much  was  used,  x^sec*  =  Xj  ,  and  form  a 

all  activities  j 


new  estimate  of  the  {K^}: 
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,New 

‘i 


,01  d 

‘i 


2x"sed 


Used  n 

X.  +  xf 


(5-13) 


If  K^ew  =  K?ld,  the  equilibrium  is  known.  If  not,  go  back 
to  Step  2. 
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6.0  SINGLE  TIME  PERIOD  OPTIMIZER 

The  single  time  period  optimizer  (subroutine  ACTVTY  in  the  model) 
optimizes  the  level  of  production  and  consumption  activities,  calculates 
the  rental  values  for  capital,  and  calculates  the  level  of  investment 
(production  of  capital).  Inputs  include  capital  inventory  available 
for  the  production  activities,  capital  value  from  the  future  (Phase  1) 
or  investment  activities  (Phase  2),  depreciation  rates  for  capital, 
translog  parameters  (see  Assumption  2  of  Sec.  2),  minimum  consumption 
activity  levels,  and  the  direct  requirements  matrix.  Outputs  include 
the  rental  values  for  capital  and  labor,  the  activity  levels,  and  the 
consumable  values.  The  logic  flow  of  subroutine  ACTVTY  is  diagrammed 
in  Fig.  6-6. 


6.1  PRODUCTION  ACTIVITIES 

Each  consumable  i  is  produced  by  a  separate  activity  using  capital 
K.j ,  labor  L^,  and  a  specified  fraction  of  each  of  the  N  consumables 
per  unit  activity.  If  more  than  one  type  of  capital  is  used,  Ki  refers 
to  the  "effective"  capital  produced  by  the  combination  (see  Sec.  4.2.4). 


The  capital  used  in  the  production  of  consumable  i  is  the  entire 
inventory  of  capital  available.  Obviously  one  cannot  use  more  capital 
than  available,  and  using  less  would  cause  its  value  to  drop  to  zero 
relative  to  labor.  Therefore,  from  the  translog  function. 


>K,Ki 


-6..  -Bi 

1  +  v  L 

V1 


■l/&i 


(6-1) 


the  level  of  production  activity  is  dependent.on  labor  alone.  But  how 
much  labor  should  be  applied?  Section  2  answered  that  question  through 
Eq.  2-10: 


-Bi/a+a,) 


(6-2) 


Li  -  Ki 


J 


where  the  parameters  were  described  in  Sec.  2.  The  important  variables 
in  Eq.  6-2  are  the  A!,  the  gross  profit  from  producing  consumable  i  which 
is  available  for  renting  the  capital  and  labor  used  to  produce  it: 

A)  5  A,  -  Z)  ojJAJ 

J 

By  combining  Eqs.  6-1  and  6-2,  Al  may  be  related  directly  to  the  activity 
level  X*?: 


Equations  6-1  through  6-3  are  the  crux  of  the  single  time  period 
optimization.  The  process  starts  as  follows: 

1.  Assume  some  level  of  labor,  L. ,  for  each  production 

activity.  This  labor  is  usually  the  labor  required  to 
produce  the  same  level  of  activity  as  achieved  for  this 
consumable  and  time  period  last  pass. 

From  calculate  x[?  from  Eq.  6-1. 


2. 


3. 


From  X?  calculate  the  ratio  A„  /{e.y^A'.)  for  each 
consumable  from  Eq.  6-3. 

Sum  each  L..  to  find  the  total  labor  used  and  calculate 
A^  from  the  resulting  consumption  of  leisure: 


r Leisure 

A  |  ~  ~  h 

L  /  V-"  .  \  Leisure 


1 


Find  A,,  from  Eq.  2-8: 

S' 


(6-4) 


Vi 


W  W 


(6-5) 


6.  From  the  ratio  found  in  step  3  and  A.,  ,  calculate 
the  A! .  1 

At  this  point  the  set  is  known,  but  the  individual  A.  are 
not.  However,  through  the  definition  of  the  A'-  in  matrix  notation, 

A'  =  (I  -  o)A 


the  (A^ }  may  be  found  by  inverting  the  matrix  (I  -  a).  From  the  A^ 
the  consumption  activities  may  be  calculated. 


6.2  CONSUMPTION  ACTIVITIES 

Once  the  A.  are  known,  the  cost  per  unit  of  consumption  activity  i 
may  be  calculated.  If  the  consumption  activity  i  uses  of  consumable  j, 
then  the  cost  per  unit  of  activity  is: 


Ci  =  ?aijAj 


(6-7) 
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The  activity  level  is  determined  by  Eq.  2-2 


i 


x 


c 

i 


(6-8) 


6.3  INVESTMENT  ACTIVITIES 

Investment  activities  are  calculated  in  one  of  two  ways  depending 
on  whether  the  program  has  just  started  or  is  trying  to  refine  its 
trajectory. 

6.3.1  Phase  One 

In  phase  one  of  the  program,  the  level  of  investment  activities 
is  found  by  comparing  the  cost  of  producing  capital  with  its  expected 
worth  to  the  future.  Assuming  the  program  has  completed  at  least  one 
pass,  there  exists  a  time  stream  of  rental  values  for  the  various  capital. 
These  may  be  propagated  backwards  in  time,  accounting  for  discount  rates, 
depreciation,  etc.,  as  described  in  Section  7.2,  to  yield  a  value  for 
capital  to  each  time  period.  The  cost  of  producing  the  capital,  which 
uses  only  consumables  in  its  production,  may  be  found  in  a  way  identical 
to  Eq.  6-7. 


The  problem  arises,  however,  that  if  the  particular  time  period 
has  a  cost  of  production  greater  or  less  than  the  worth  of  the  capital, 
it  will  want  to  produce  zero  or  infinite  amounts  of  the  capital.  This 
problem  is  essentially  solved  by  assuming  the  value  of  the  capital  is 
correct  only  if  the  given  time  period  makes  capital  in  reasonable  quan¬ 
tities.  Less  production  increases  the  future  value  of  the  capital,  and 
more  production  reduces  the  future  value  of  the  capital.  The  functional 
nature  of  this  process  is  described  in  terms  of  two  parameters,  RTGT 
and  Rpy-p  Rpyj  is  the  capital  resource  level  one  gestation  time  in  the 
future  if  nothing  is  produced.  RjQj  is  the  target  resource  level,  i.e., 
it  is  the  level  of  inventory  which  existed  one  gestation  time  in  the 
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future  during  the  last  pass.  The  difference,  -  Rpyj.  is  what  should 
be  produced  this  time  period  assuming  cost  equals  value. 1 

Once  Rpyy  and  R^y  are  known,  xj  is  found  through  the  equation; 

X1  -  RTGT  '  RFUT  *  £i  1n  ^Ai/Cj)  ,  . 

i  At  v  ' 

where  is  the  value  of  the  capital  propagated  from  the  future,  and 
s.  is  a  scale  factor  essentially  equal  to  the  equilibrium  activity  level 
x|eq.  The  effect  of  Eq.  6-9  is  to  produce  negative  feedback.  If  the 
investment  activity  is  too  small,  rental  values  will  be  high  in  the 
future  causing  A.  to  be  large  next  pass,  causing  the  time  period  to 
produce  more  than  it  did  last  pass.  This  in  turn  causes  the  inventory 
in  the  future  to  rise  and  rental  values  to  fall,  lowering  a|  for  the 
next  pass. 

While  this  process  sounds  as  if  it  should  cause  swift  convergence, 
it  in  fact  causes  too  much  feedback,  resulting  in  an  oscillatory  state 
wherein  every  other  pass  results  in  too  little  investment  followed  by 
too  much  investment.  This  effect  is  cured  by  averaging  inventories 
from  pass  to  pass  until  a  pseudo-stable  trajectory  arises,  and  then 
proceeding  into  phase  two. 

For  the  first  pass  of  the  program  RTGT  1s  defined  as  the  equili¬ 
brium  resource  level  while  A*  is  the  equilibrium  capital  value.  These 
values  cause  the  model  to  invest  to  such  a  degree  that,  while  not 
producing  an  optimal  trajectory,  it  at  least  produces  one  which  brings 
the  economy  to  equilibrium. 

1 In  reality,  more  needs  to  be  produced  because  of  depreciation.  There¬ 
fore  the  model  uses  Rpyy  2  RqUT(1  -  At/F)  +  Rp^At/F) , 

where  F  =  e*Px  (1  -  e'(P+d^AtV(P+d)  (See  Sec.  7.1). 
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6.3.2  Phase  Two 

In  phase  two  the  investment  trajectory  is  approximately  known. 

At  this  stage  the  single  time  period  optimizer  accepts  the  investment 
activities  as  given,  regardless  of  the  cost  of  production  or  the  future 
value  of  the  capital.  These  activities  are  adjusted  up  or  down  external 
to  the  single  time  period  optimizer  through  the  use  of  an  Everett  type 
algorithm. 

This  algorithm  operates  in  the  following  way.  Consider  a  single 
activity  X*.  On  the  previous  pass  the  activity  had  a  cost  of  production 
equal  to  cJ.  After  the  rental  values  have  been  propagated  back  in  time 

1  T 

a  value  A..  results.  For  the  next  pass  X..  is  reset  by: 


X 

►— t 

II 

X 

— *•  * — t 

»— » 

+ 

ro 

C1 

Li 

<Ai 

(6-10) 

or 

*i  -*}u  -  o. 

Li 

>V 

The  parameter  e  is  reset  each  time  by 

e  =  e(l  +  SF) ,  for  S  =  +1 

or  e  -  e(l  -  6),  for's  -  -1  (6-11) 

I  I 

where  S  is  defined  as  the  direction  of  change,  i.e,  S  =  +1  if  C.  >  A. 

II  11 

both  this  pass  and  last  pass  or  C.  <  A.  both  this  pass  and  last  pass. 

It  '  *  I  I 

S  =  -1  if  C.  >  A.  for  one  of  the  two  passes  and  Ci  <  A.  for  the  other. 

The  parameters  <5  and  F  are  constants.  In  the  model  they  are  set  to  0.6 

and  0.5  respectively.  This  algorithm  very  quickly  and  efficiently 

narrows  the  investment  trajectory  to  the  optimum.  Inventories  are  not 

averaged  during  phase  two. 
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6.4  TRADE 

The  current  ACDA  version  of  DYNEVAL  does  not  model  trade  explicitly. 
Some  mechanism  is  needed,  however,  to  account  for  the  "undifferentiated 
imports"  required  in  some  production  and  investment  activities.  This 
is  handled  by  summing  all  activities  to  determine  the  amount  of  undiffer¬ 
entiated  imports  required,  and  setting  the  trade  activity  at  the  proper 
level.  All  other  commodities  are  imported  or  exported  at  this  level 
according  to  the  proportions  in  which  they  are  imported  or  exported  in 
the  current  economy. 

6.5  CONVERGENCE 

Once  the  production,  consumption,  trade  and  investment  activities 
are  known,  a  comparison  may  be  made  between  the  amount  of  a  given  consum¬ 
able  which  is  produced  and  the  amount  which  is  used  or  consumed.  These 
two  values  should  be  equal.  In  general,  however,  they  will  not  be,  and 
a  scheme  must  be  developed  to  bring  them  together.  This  is  accomplished 
by  using  once  again  an  Everett  type  of  algorithm  on  the  A!  *  of  each 
consumable.  When  production  exceeds  consumption  A!  is  reduced  and  when 
consumption  exceed  production  A!  is  increased. 

When  A*  is  changed,  all  other  variables  will  also  change.  One 
parameter,  however,  is  less  sensitive  than  most.  This  is  the  rental  value 
of  labor.  Using  Eq.  2-10b  with  the  new  value  of  Al  and  the  old  value  of 
\L  provides  a  new  estimate  for  L^.  One  may  then. proceed  as  discussed 
in  6.1  above. 


llt  is  important  that  the  Everett  algorithm  be  applied  to  the  Aj  and  not 
the  individual  A,  Doing  the  latter  results  in  wild  swings  in  the  pro¬ 
duction  activities  since  Ai  can  easily  become  negative  resulting  in 
zero  production.  A  higherA^  must  result  in  increased  production. 

A  higher  A-j  need  not  do  so  if  the  other  increase  as  well. 


This  process  usually  converges  quickly  and  efficiently.  When 
there  has  been  a  large  destruction  of  capital,  however,  the  early  time 
periods  may  have  trouble  converging.  This  occurs  because  one  or  more 
industries  is  working  near  capacity  and  yet  its  keeps  getting 
signals  to  increase.  The  industry  cannot  increase  its  output  so 
costs,  including  labor  costs,  increase  until  consumers  cannot  purchase 
anything  and  the  system  reverses  itself.  To  prevent  the  model  from 
spending  too  much  of  its  time  on  a  problem  which  will  be  changed  anyway 
after  the  future  sends  its  information  back  through  time,  a  limit  of 
20  iterations  is  imposed  on  the  single  time  period  optimizer. 
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7.0  PROPAGATION  THROUGH  TIME  OF  CAPITAL  AND  VALUE 

Figure  7-1  illustrates  the  temporal  structure  of  the  model.  Each 
time  period  K  is  of  equal  length,  At.  The  array,  T(K),  stores  the  time 
at  the  middle  of  time  period  K,  [i.e.,T(K)  =  (K  -  l)At  +  At/2].  Tq(K)  is 
the  time  at  the  start  of  time  period  K.  The  activities,  X,  and  the 
shadow  (rental)  values,  X,  are  represented  as  constant  rates  within  time 
period  K: 


Resource  inventories  are  considered  continuous  in  time,  with  the  values 
at  T(K)  stored  in  array  R(K):1 


*In  this  section  capital  resources  will  be  designated  as  R  to  conform 
to  the  model  representation  and  to  avoid  confusion  with  the  time  index  K. 
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Figure  7-1.  Temporal  Structure  of  the  Model 


The  shadow  value  A(K  +  1)  represents  the  worth  of  the  resource  which 
time  period  K  passes  to  the  future  measured  in  present  day  value.  It 
does  not  represent  the  worth  to  time  period  K  +  1.  Rather,  A  is  used 
by  time  period  K  in  determining  the  amount  it  should  produce,  and  so 
reflects  the  discounted  future  stream  of  X's  encountered  after  the 
gestation  time  t. 

The  resource  level  at  the  center  of  time  period  K,  R^(K),  is  the 
inventory  available  to  time  period  K.  Strictly  speaking,  an  integrated 
average  of  Ri  would  be  a  better  approximation.  However,  for  gestation 
times  t  less  than  At  the  averaging  would  be  incorrect  anyway.  For 
example,  if  t  =  J^At,  the  following  picture  arises: 

„  ^  ~  Future  course  of 


x '  inventory  level 

Decay  if  no  activity  in 

^  ! 

i 

_  time  period  K 

l 

i 

•  i  t 

-2  K-l 

T(K) 

K  K+l 

,  „ 

K+2 

At  the  beginning  of  time  period  K,  the  activity  level  X(K)  has  not  been 
determined,  and  the  resource  level  at  the  end  of  time  period  K  reflects 
the  natural  depreciation  of  the  resource.  In  this  example,  using  R(K) 
is  obviously  a  better  estimate  of  the  available  resources  than  an  actual 
integration  of  the  curve  would  yield. 

7.1  RESOURCE  PROPAGATION 

After  the  X*  *  are  determined,  the  resources  produced  are  propagated 
through  time.  This  occurs  in  subroutine  TIMPER.  As  resources  propagate, 

^X j ( K)  here  refers  to  the  rate  of  production  of  resource  R.  during  the 
Kth  time  period  in  units  of  resource/unit  time. 
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they  depreciate  with  a  natural  depreciation  rate  d.  Furthermore,  resources 
are  normalized  to  population  so  as  the  population  grows  at  a  rate  P, 
resource  inventories  are  reduced.  As  the  resource  propagates,  each  value 
of  R(K‘  >  K)  is  increased  by  an  amount  AR(K'),  where  AR(K')  is  determined 
in  the  following  way. 

Resources  begin  to  come  on  line  at  time  T$tart  =  T0(K)  +  t 

*  T(K)  -  +  T.  For  time  t  _  T$tart»  AR  =  0.  No  resource  is  yet 

available  due  to  the  gestation  time  t.  Therefore,  for  K'  such  that 
T(K' )  _<  Ijtart’ 

AR(K')=0  (7-D 

Resources  finish  coming  on  line  at  time  T$top  =  TQ(K  +  1)  +  t 
=  T(K  +1)  -  +  t.  For  time  t  between  Tstart  and  TstQ- ,  AR  reflects 

the  incoming  resources  which  have  come  on  line  by  time  t.  Therefore, 
for  K'  such  that 


Tstart  <  W)  1  T$top, 


I  -Px  fT(<  )_Tstart  e-(p+dK  dt 
AR(K')  =  Xl(K’)e  1  e 


f 


.  x‘(K')  e"Px 


-(P+d)[T(K')-Tstart] 


P  +  d 


(7-2) 


where  d  =  depreciation  rate 

P  =  population  growth  rate 


*xj(K)  here  refers  to  the  rate  of  production  of  resource  R .  during  the 
Kth  time  period  in  units  of  resource/unit  time. 
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For  time  t  1  Tstop,  ar  reflects  the  depreciated  value  uf  all  the  resources 
produced  in  time  period  K: 


,  x  ,  v  .p  i  .  e-(P+d)At  -(P+d)[T(K')-T  ] 

AR(K  )  =  X( K)  e  P  ± ^ -  e  stop  (7-3) 

/  t  \ 

depreciation 


population 
growth  during 
gestation 


P  +  d 

t 

depreciation  between 
time  coming  available 
and  time  when  all 
units  come  on  1  ine 


after  all  units 
are  on  1 ine 


Equations  7-1  through  7-3  may  be  combined  as  follows: 


Note  that  for  x  <  ,  AR(K)  f  0.  This  is  allowed  for  all  time  periods 

except  K  =  1,  since  by  definition  R(l)  is  the  starting  resource  and  if  it  is 
allowed  to  increase,  the  next  pass  of  the  model  will  in  effect  be  solving 
a  different  problem  with  a  higher  starting  resource.  Therefore,  the  length 
of  each  time  period  must  be  such  that  At  £  2x.  The  smaller  At  is,  the 
longer  each  pass  of  the  model  takes  since  more  time  periods  are  required, 
but,  in  general  fewer  passes  are  required  due  to  increased  accuracy. 

i 

7.2  VALUE  PROPAGATION 

This  section  describes  how  the  A(K  +  1)  are  derived  for  the  next 
pass.  After  a  complete  pass,  subroutine  LAMBAC  propagates  the  rental 
values  X  back  in  time  to  define  a  future  worth,  A(K  +  1),  to  time 
period  K  for  its  capital.  Since  it  is  assumed  that  the  activity  level 
X*(K)  is  constant  throughout  the  time  period  and  that  the  decision 
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maker  decides  this  level  at  the  beginning  of  the  time  period,  the  basic 
equation  for  A(K+1)  is: 


A(K  +  lJX^KjAt'  =  discounted  rents  on  amount  produced  (7-5) 

future 

The  value  At'  represents  the  fact  that  the  value  A(K  +  1)  is  also 
discounted  through  the  time  period,  i.e., 

.  »•»  - 

P 

where  p  is  the  discount  rate.  Throughout  this  discussion  it  will  be 
assumed  that  Fig.  7-1  applies.  The  index  KK  in  Fig.  7-1  refers  to  the 
first  index  after  K  such  that  gestation  is  completed: 


KK  =  INT 


M  +..H-+  t  +  j  5 
At  1,5 


In  the  (KK-2)th  time  period,  resources  start  to  come  on  line 
and  continue  (in  this  time  period)  for  a  time  Ax  (see  Fig.  7-1).  There 
fore,  the  first  component  of  A(K  +  1)  is: 


A*(K  +  UX^KjAt'  =  A(KK  -  2)e 


renta 


value 


-(P*ph  /?„.(„  1  -f}™*  e'0t  dt 

♦  Jo  | 


P  +  d 


discount  and 
population  growth 
during  gestation 


amount  available 
between  x>t>x+Ax 


discount 

beyond 

gestation  time 


=  A(KK  -  2)  e"(p+p)T  x!(K) 


e-(P+d+p)Ax  _  1  e-pAx  _  1 


P  +  d  +  p 


(7-6) 


In  the  (KK-l)th  time  period,  resources  continue  to  come  on  line 
for  a  time  at  At  -  At.  They  all  then  depreciate  for  the  rest  of  the 
time  period  Ax.  Thus: 

A2(K  +  1 )X! (K)At '  =  A(KK  -  I)  e'(P+p)T  X!(K) 


/•A  t 

!  .  e_(P+d)t 

'Ax 

P  +  d 

e“pt  dt 


l  _  e-(P+d)At 


P  +  d 


•pAt  . 


£ 

l 


e-(P+d+p)t  dt 


-  A(KK  -  1)  e~(P*p)l  X1  (K) 
P  +  d 


( P+d+p )At  _  g-( P+d+p) Ax 


P  +  d  +  p 


e“pAx  -  e'pAt  e~pAt  -  e_(p+d+p)At  x  .  e-(P+d+p)Ax 
p  +  P  +  d  +  p 

For  the  rest  of  time,  the  resources  continue  to  depreciate: 

A3(K  +  1  )X 1  ( K)At *  =  e  ”^P+P^TXI(K)  X(KK  +  k)  e'pAt  L  ~pe^d 


(7-7) 


-(P+d+p)  (Ax+kA  ) 


1  -  e 


•(P+d+p) At 


P  +  d  +  p 


(7-8) 
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Similar  expressions  may  be  derived  for  A(K  +  2).  If  this  is  done 
it  is  seen  that  A ( K  +  1)  may  be  related  to  A ( K  +  2)  and  the  three  rental 
values,  X(KK  -  2),  X(KK  -  1),  and  A ( KK ) : 

A ( K  +  1)  =  A (K  +  2)  e-(p+d+p)At 


e-(P+p)x 
+  (P"+  d )  At 1 

+  \(KK  -  1) 


A  (KK  -  2) 

e~pAT  n  +  e 


.  e-pAT 

P 

-(P+d+p)  At 

p 


1  .  e-(P+d+P)Ax 
P  +  d  +  p 


1  .  e-(p+d+p)At 


e-pAt;^  _  g- (P+d+p)  At -j  +  j-g-( P+d+p)  A  t  _  e-(P+d+p)Ax-| 

P  +  d  +  p 


+  A ( KK) 


e-pAt  j-e-( P+d+p )At  _  g-(P+d+p)At-| 
P  +  d  +  p 


+  e-( P+d+p) At 


(e-pAt  _  e-pAx) 
p 


(7-9) 


Equation  7-9  is  used  by  the  mo.del  in  subroutine  LAMBAC  to  derive 
the  shadow  values  A(K  +  1)  used  in  the  next  pass.  Note  that  at 
equilibrium  A(K  +  1)  =  A(K  +  2)  and  X(KK  -  2)  =  X(KK  -  1)  =  X(KK). 

In  this  case,  Eq.  7-9  reduces  to: 


x  e-(P+d+p)^ 
P  +  d  +  p 


(7-10) 


which  is  the  equilibrium  value  used  in  Sec.  5.1. 
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APPENDIX  A 
USER'S  MANUAL 
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A.  1  INTRODUCTION 

This  appendix  explains  how  to  use  the  Dynamic  Economic  Values  Model 
(DYNEVAL).  It  is  assumed  that  a  data  file  exists  containing  the  primary 
data  describing  the  current  economy.  This  file  is  most  easily  created 
through  the  use  of  the  data  aggregator  AGGRAT  which  aggregates  DSA 
format  economic  data  into  economic  sectors  specified  by  the  user  and 
puts  the  aggregated  data  onto  a  file  in  a  form  accepted  by  the  model. 
AGGRAT  is  described  in  Appendix  B. 

Besides  the  data  describing  the  current  economy,  which  includes 
depreciation  rates,  gestation  times  for  capital  formation,  etc.,  spe¬ 
cific  variables  describing  the  disruption  or  change  to  the  economy 
are  needed.  These  variables  are  unique  to  each  study  and  need  to  be 
entered  when  the  model  is  run.  The  following  is  a  step  by  step 
description  of  how  these  variables  are  entered,  what  options  are 
available  to  the  user,  and  what  the  final  output  of  the  model  (termi¬ 
nal  and  line  printer)  looks  like.  A  sample  listing  of  line  printer 
output  is  given  in  Sec.  A. 5. 

A. 2  START-UP 

Upon  initiating  execution  of  DYNEVAL,  the  following  question  is 

asked: 


IS  THIS  A  RESTART? 

The  model  has  a  restart  capability  which  allows  the  user  to  stop  the 
model  and  restart  it  where  it  left  off  with  changes  in  print  and  plot 
switches  (see  Sec.  A. 3).  The  restart  capability  is  especially  useful 
when  peculiarities  in  the  summary  outputs  aresnoticed  and  the  user 
wants  to  get  detailed  prints  or  plots  of  the  next  couple  of  passes. 
This  question,  like  all  yes/no  questions,  may  be  answered  with  YES(NO) 
or  simply  Y(N). 
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OYNEVAL  next  asks: 


FULL  PRINT  OF  EACH  PASS  (Y)  OR 
SUMMARY  ONLY  (N) 

This  question  refers  to  line  printer  output.  The  model  calculates  the 
optimal  economic  trajectory  through  an  iterative  technique.  A  summary 
of  each  pass  is  printed  on  both  the  terminal  and  the  line  printer  (see 
Sec.  D).  If  the  user  answers  this  question  with  YES,  however,  a  time 
period  by  time  period  listing  of  important  parameters  is  also  printed 
on  the  line  printer.  A  description  of  this  output  is  given  in  Sec.  A. 4 

Besides  numerical  output,  the  model  is  also  capable  of  graphical 
summaries.  The  next  question  asked  by  DYNEVAL  is: 

IS  OUTPUT  FOR  PLOTS  DESIRED?  (YES,  NO,  OR  END). 

A  YES  or  END  answer  opens  (creates)  a  plotting  file  and  sets  the  plot 
switch.  If  the  user  wants  to  see  how  selected  variables  change  from 
pass  to  pass  or  wants  to  make  graphs  at  various  stages  of  convergence, 
he  should  answer  with  YES.  If  the  user  does  not  care  about  the  con¬ 
vergence  process  but  would  like  time  period  plots  of  the  final  economic 
trajectories,  he  should  answer  END.  It  should  be  noted  that  if  the 
plot  option  is  selected  the  model  merely  outputs  relevant  data  onto  a 
file.  A  separate  program  is  required  to  convert-  the  data  into  plots. 
Such  a  program  (PLOTER)  has  been  written  using  the  Tektronix  software, 
PLOT-10.  Should  a  different  software  package  (e.g.,  DEC'S  DSSPLA)  be 
desired,  a  new  program  will  have  to  be  designed. 

Once  the  output  of  the  model  has  been  determined,  DYNEVAL  queries 
the  user  concerning  the  parameters  for  the  study.  The  first  data 
asked  for  is: 

INPUT  PASS  LIMIT  AND  STOPPING  CRITERION 
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The  model  will  continue  its  operation  until  either  convergence  is 
achieved  or  it  has  completed  a  maximum  number  of  passes.  The  pass 
1 imit  is  the  maximum  number  of  passes  the  user  wants  the  model  to 
complete.  A  value  of  50-100  is  normally  sufficient.  The  stopping 
criterion  is  the  value  of  F0M2,  the  second  f igure-of-merit  printed  in 
the  summary  of  each  pass,  which  indicates  that  convergence  is  satis¬ 
factory.  The  various  f igures-of-meri t  are  discussed  in  Sec.  A. 3.  A 
value  of  0.01-0.05  for  the  stopping  criterion  is  usually  sufficient. 

If  this  is  a  restart,  the  program  proceeds  where  it  left  off. 
Otherwise,  DYNEVAL  next  asks  for  the  extent  of  the  run: 

INPUT  NUMBER  OF  TIME  PERIODS  DESIRED 

The  number  of  time  periods  times  the  length  of  each  time  period  should 
be  large  enough  so  that  the  economy  has  time  to  reach  its  equilibrium 
state.  Forty  to  fifty  years  is  sufficient  for  most  studies.  After 
the  number  of  time  periods  has  been  input,  the  model  is  able  to  cal¬ 
culate  the  total  storage  requirement  for  the  study.  It  informs  the 
user  how  much  storage  he  is  using  and  how  much  is  available: 

STORAGE  USED  =  8732,  STORAGE  AVAILABLE  =  20,000 

If  too  much  storage  is  required,  the  following  message  will  be  printed 
on  the  terminal : 

INCREASE  DIMENSION  OF  ARRAY  W  IN  MAIN  PROGRAM  TO  AT 
LEAST  (amount  required)  --  PROGRAM  ABORTED 

The  user  will  have  to  increase  the  dimension  of  array  W  and  change  the 
value  of  parameter  NW  in  the  main  program  and  recompile.  Ample  storage 
is  allowed  for  most  studies. 
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Assuming  the  storage  requirement  is  smaller  than  the  available 
storage,  DYNEVAL  next  asks  for  the  length  of  each  time  period: 

INPUT  LENGTH  OF  TIME  PERIOD  (YEARS) 

For  a  given  length  of  time  over  which  the  model  optimizes,  greater 
accuracy  with  fewer  passes  is  achieved  with  many  time  periods  of  short 
duration.  More  data  storage  and  longer  time  per  pass  is  also  required, 
however,  so  the  user  needs  to  trade  off  the  costs  and  benefits  of  the 
division.  In  no  case  may  the  time  period  length  be  greater  than  twice 
the  smallest  gestation  time.  If  this  occurs,  DYNEVAL  will  warn  the  user: 

TIME  PERIOD  LENGTH  MAY  NOT  BE  GREATER  THAN  (NO,  OF  YEARS), 

I.E.,  2.*  MINIMUM  GESTATION  TIME  OF  (SMALLEST  GESTATION 
TIME)  INPUT  LENGTH  OF  TIME  PERIOD  (YEARS) 

The  model  next  asks  the  user  about  the  basic  parameters  of  the 
study.  It  remembers  the  parameters  used  for  the  last  study,  but  the 
user  may  change  any  he  wishes.  If  the  user  would  like  to  see  the  current 
values  of  the  basic  parameters  he  may  either  answer  YES  to  the  following: 

WOULD  YOU  LIKE  BASIC  PARAMETER  ARRAYS 
PRINTED  ON  LP?  (Y  or  N), 

which  will  output  the  basic  arrays  to  the  line  printer,  or,  if  this 
question  is  answered  with  NO,  he  may  have  specific  arrays  displayed  on 
the  terminal  by  answering  YES  to: 

WOULD  YOU  LIKE  ANY  OF  THE  BASIC  PARAMETER  ARRAYS 
DISPLAYED?  (Y  or  N) 

Assuming  the  user  answers  YES,  DYNEVAL  answers  with: 

INPUT  ARRAY  NAME 

TYPE  'M'  TO  DISPLAY  MENU 
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Entering  an  M  on  the  terminal  shows  the  user  what  the  "basic  parameter 
arrays"  are  and  how  to  access  them. 

MENU:  Z  =  MINIMUM  CONSUMPTION  ACTIVITY 

B  =  UTILITY  EXPONENT  FOR  CONSUMPTION  ACTIVITY 

G  =  RATIO  OF  POST-ATTACK  UTILITY  COEFFICIENTS  TO  TODAY'S 
UTILITY  COEFFICIENTS 

BETA  =  TRANSLOG  EXPONENT  FOR  CAPITAL/LABOR  TRADE-OFFS 
(IF  .LT.  -1.0,  BETA  IS  CALCULATED  TO  GIVE  THE 
LIMITING  PRODUCTIVITY  A  VALUE  OF  ABS  (INPUT 
VALUE)  (See  Sec.  4.2.3} 

BET2  =  TRANSLOG  EXPONENT  FOR  CAPITAL/ CAPITAL  TRADE-OFFS 
(ONLY  FOR  THOSE  PRODUCTION  ACTIVITIES  USING 
2  TYPES  OF  CAPITAL) 

F  =  FRACTION  OF  CAPITAL  SURVIVING 

M  =  REDISPLAY  MENU 

E  =  END  OF  REQUESTS 

INPUT  ARRAY  NAME 

TYPE  'M'  TO  DISPLAY  MENU 


After  perusing  the  data,  the  user  may  wish  to  change  some  of  it: 

WOULD  YOU  LIKE  TO  CHANGE  ANY  OF  THE  'DATA  (Y  or  N) 

If  the  user  answers  with  YES,  he  may  change  any  or  all  of  the  data 
provided  in  the  above  menu: 


INPUT  ARRAY  NAME,  DATA  CHANGE,  AND  INDUSTRY  NUMBER 
(OR  RANGE  OF  INDUSTRY  NUMBERS) 

TYPE  'M1  TO  DISPLAY  MENU 


. 
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The  data  arrays  are  dimensioned  by  the  number  of  consumables,  number  of 
capital  items,  number  of  consumption  activities,  etc.  If  the  user  tries 
to  change  a  piece  of  data  outside  its  normal  range,  DYNEVAL  will  inform 
the  user  of  his  error: 

EXCUSE  ME? 

RANGE  MUST  BE  BETWEEN  1  AND  5 
(values  for  example  only) 

When  all  changes  are  complete,  the  user  has  the  option  of  once  again 
seeing  the  data  either  printed  on  the  line  printer  or  displayed  on  the 
terminal.  If  all  is  in  order,  the  run  begins  with  the  calculation  of 
the  final  equilibrium  state.  When  this  is  completed,  the  model  tells 
the  user  how  long  the  calculation  took  by  printing: 

CPU  TIME  TO  CALCULATE  EQUILIBRIUM  =  (NUMBER)  CPU  SECONDS 
FOR  (NUMBER)  ITERATIONS 

on  the  terminal.  Following  this  calculation,  the  optimal  trajectory 
to  this  equilibrium  state  is  calculated  through  the  iterative  technique 
of  Dynamic  Lagrange  Programming. 

A. 3  RESTART 

If,  for  any  reason,  the  user  wishes  to  stop  the  run,  he  need  only 
type  any  character  on  the  terminal.  At  the  beginning  of  each  pass,  the 
program  checks  to  see  if  the  terminal  input  buffer  is  empty.  If  it 
is  not,  the  model  responds  with: 

EXCUSE  ME? 

YOU  HAVE  INTERRUPTED  ME 

TYPE  "C"  TO  CONTINUE,  "S"  TO  STOP 
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The  user  may  then  stop  the  program  by  typing  on  "S".  All  files  will 
be  closed  and  the  program  will  end. 

The  user  may  then  check  the  output  by  printing  file  PRINT.DAT.  If 
he  wishes  to  continue  the  run  later,  he  may  simply  reexecute  EREC  and 
answer  YES  to  the  question 

IS  THIS  A  RESTART? 

(see  Sec.  B).  All  printer  output  in  a  restart  is  appended  to  the  current 
version  of  PRINT.DAT  so  the  initial  output  is  not  lost. 

In  normal  operation,  a  restart  is  also  possible  if  the  program 
fails  for  some  reason,  e.g.,  because  of  a  system  crash,  as  long  as  the 
failure  did  not  occur  while  the  insurance  file  was  being  written.  The 
insurance  file  is  rewritten  after  each  pass  and  the  output  file  (PRINT.DAT) 
is  closed  to  save  the  output  and  then  reopened  in  the  append  mode.  To 
be  safe,  however,  it  is  advisable  to  stop  the  program  with  a  FORTRAN 
interrupt  as  described  above. 

A. 4  TERMINAL  OUTPUT 

Throughout  the  course  of  the  run,  a  summary  of  each  pass  is  dis¬ 
played  on  the  terminal,  e.g.: 

PASS  5;  24.320  CPU  SEC;  508  LOOPS;  FOM  =  0.967E-02  0.179  1.000 

The  above  summary  is  for  the  fifth  pass  which  took  24.32  cpu  sec  to 
complete.  The  number  of  loops  refers  to  the  cumulative  total  of  single 
time  period  optimization  iterations.  There  is  a  maximum  number  of 
iterations  for  any  given  time  period  (nominally  20).  The  three  figure- 
of-merit's  are  defined  for  each  pass  as: 
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FOM(l)  = 


where 


V**ik>2 


F0M(2)  = 


FOM(3)  = 


k=l  i=l 


N  M 

E  EAik  xtk 

k=l  i=l 


(la) 


(lb) 


(lc) 


N  =  Number  of  time  periods 

M  =  Number  of  capital  resources 

Tf^  =  Average  inventory  of  capital  i, 

time  period  K.  *  i(R.|k-  +  R.jJlast  pass)) 

AR^  =  Percentage  difference  in  Rik 

(ARik  =  (R-k  -  R1k(last  pass))/R.k) 

AA^k  =  Percentage  difference  between  value  of  capital 

and  effective  cost  of  producing  capital 

/ 

Aik  =  Effective  cost  of  producing  capital 

Xik  =  Amount  of  capital  i  produced  in  time  period  k 

ei<k  =  Everett  parameter  (phase  two  only). 
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The  first  f igure-of-merit  is  a  measure  of  how  closely  the  capital  inven¬ 
tory  trajectory  is  being  reproduced  each  pass.  It  is  the  square  root 
of  the  average  square  percentage  change  in  the  inventory,  weighted  by 
the  average  size  of  that  inventory.  The  second  f igure-of-merit  shows 
how  closely  the  value  of  the  capital,  found  by  propagating  the  rental 
values  back  in  time,  matches  the  cost  of  producing  the  capital.  It  is 
simply  the  RMS  value  of  the  percentage  difference  of  the  two  weighted 
by  the  size  of  the  investment  activity.  The  third  figure-of-merit 
is  a  measure  of  how  well  the  Everett  algorithm  of  phase  two  (see 
Sec.  6.3.2)  is  performing.  The  change  in  investment  activity  level, 
e^,  is  weighted  by  the  total  cost  of  capital  produced. 

Phase  two  of  the  model  is  started  whenever  one  of  the  following 
is  achieved: 


1.  FOM  (1)  <  .005 

2.  K  >  20 

3.  K  >  N/2 

where  K  is  the  current  pass  number.  When  phase  two  begins,  the  following 
message  is  printed  on  the  terminal: 

***BEGINNING  PHASE  TWO*** 

The  user  will  notice  that  after  each  summary  print  on  the  terminal 
the  cursor  will  remain  at  the  end  of  the  line  for  a  short  period  and 
then  return  to  the  beginning  of  the  next  line.  The  cursor  return  is 
a  signal  that  the  insurance  file  (INS. DAT)  has  been  written.  This 
insurance  file  is  used  by  the  model  whenever  the  user  specifies  that 
the  run  is  a  restart.  It  is  rewritten  at  the  end  of  each  pass. 

A. 5  LINE  PRINTER  OUTPUT 

A  sample  of  the  line  printer  output  is  given  in  Sec.  A. 5.1  below. 
The  output  from  the  line  printer  begins  by  describing  the  basic  parameter 
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arrays  for  this  run.  These  consist  of  the  length  of  the  run,  number  of 
consumables,  etc.,  minimum  activity  levels,  utility  function  elasticity 
parameters  (b),  the  ratio  of  post-attack  to  preattack  utility  function 
coefficients  (F),  the  input  elasticity  parameters  for  the  translog 
substitution  function  (8)  (see  Sec.  A. 2),  and  the  fraction  of  capital 
surviving. 

Following  the  above  is  a  listing  of  the  basic  data  arrays  for  the 
economy  under  study:  current  activity  levels,  resource  inventories, 
gestation  times,  depreciation,  growth  factors,  discounts,  direct  require¬ 
ments  matrix  for  consumables ,  the  current  capital  and  labor  usage  for 
each  production  activity,  and  a  summary  of  the  level  of  the  current 
economy  (e.g.,GNP). 

After  these  basic  arrays,  derived  current  day  values  for  functional 
parameters  are  displayed.  These  include  the  actual  elasticity  parameters 
(8),  the  utility  function  coefficients  (F),  the  rental  values  cor  capital 
and  labor  (X),  the  translog  parameters  for  the  mix  of  capital  and  labor 
(y)  and,  if  two  types  of  capital  are  used  in  a  production  activity, 
the  translog  parameters  for  the  effective  capital,  and  the  translog 
normalization  coefficients  (e). 

The  equilibrium  values  of  selected  arrays  are  displayed  next.  These 
include  the  new  utility  function  coefficients  (r),  resource  inventories, 
activity  levels,  values,  and  capital  rental  values.  A  listing  of  the 
productivity  of  each  production  industry  as  a  function  of  the  labor/ 
capital  ratio  is  also  given  as  an  aid  to  determining  those  industries 
which  are  capital  intensive  and  those  which  are  laDor  intensive.  Following 
these  arrays  is  again  a  summary  of  equilibrium  economic  levels. 

The  output  next  shows  the  level  of  capital  inventories  at  time 
zero  and  proceeds  to  duplicate  the  terminal  output  (Sec.  A. 4).  If  the 
user  specified  during  the  start-up  phase  that  he  wanted  detailed  prints, 
then  the  output  will  also  contain  a  listing  of  the  inputs  and  outputs 
as  shown  in  Fig.  A-i . 
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Figure  A-l.  Sample  Listing  of  the  Inputs  and  Outputs 


R  refers  to  the  capital  resource  inventory  (including  labor).  In 
this  case  there  are  eight  capital  items,  each  used  in  one  of  eight 
production  industries.  FLAM  refers  to  the  value  of  capital  (propagated 
from  the  future  or,  in  this  case  for  the  first  past,  set  equal  to  the 
equilibrium  value).  The  last  eight  values  of  FLAM  are  the  values  of 
the  eight  consumables.  These  values  have  meaning  only  for  the  output 
FLAM'S.  RTRGT,  RFUT,  and  XTRGT  refer  to  the  target  resource  inventory, 
i.e.,  the  inventory  which  will  occur  one  gestation  time  in  the  future 
if  no  investment  is  made,  and  the  corresponding  target  activity  level. 
These  parameters  have  meaning  and  are  printed  for  phase  one  only.  DLAMB 
is  the  rental  value  for  capital  (and  labor).  X  refers  to  the  activity 
level  of  all  activities.  EFFCST  is  the  cost  of  production  for  the 
capital.  LABOR  is  the  amount  of  labor  used  in  producing  each  commodity. 
LAMPRM  (A1)  refers  to  the  gross  profit  for  each  commodity  which  may  be 
divided  between  capital  and  labor.  Finally,  TCONS  is  the  total  amount 
of  each  commodity  used  or  consumed.  It  should  be  compared  with  the 
total  amount  X  produced  (elements  9  -  16  of  the  X  array). 

Once  the  model  converges  (or  reaches  the  pass  limit),  the  final 
results  are  printed.  These  include  the  inventory,  activity  level, 
shadow  value  and  capital  rental  value  trajectories.  Cost  of  capital 
production,  percentage  difference  with  capital  value,  and  Everett 
parameters  (phase  two)  are  also  printed  so  the  extent  of  convergence 
may  be  determined.  A  summary  of  the  economy  (gross  consumption,  in¬ 
vestment,  GNP,  etc.)  ends  the  output. 

A. 5.1  Sample  Listing  of  Line  Printer  Output 
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DATA  AGGREGATOR- 


-AGGRAT 
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The  user  will  usually  want  to  aggregate  the  economic  data  into 
sectors  of  interest.  Perhaps  he  would  like  to  keep  the  energy  sectors 
separate,  but  combine  many  manufacturing  sectors  into  one  sector.  This 
can  most  easily  be  accomplished  through  the  use  of  the  DSA  data  aggregator, 
AG6RAT,  which  also  performs  the  task  of  calculating  or  asking  for  missing 
data  such  as  depreciation  rates,  discount  rates,  etc.,  and  putting  the 
final  data  into  model  format. 

B.l  DSA  FORMAT 

A  primary  data  file  must  exist  for  AGGRAT  to  use.  This  file  consists 
of  a  series  of  records,  where  each  record  contains  four  entries  labeled 
I,  J,  K,  and  DATA.  These  four  entries  constitute  the  raw  economic  data 
in  DSA  format,  as  shown  in  Fig.  B-l.  The  index  K  refers  to  "pages"  1-5. 
Indices  I  and  J  refer  usually  to  industry  number,  and  are  the  column  and 
row  indices  for  each  page.  The  file  must  be  sorted  by  K,  J,  and  I. 

B.1.1  PAGE  1  -  PRODUCTION 

The  first  page  of  the  file  is  the  production  transactions  matrix. 

If  there  are  N  production  industries,  this  page  consists  of  N  columns 
and  N  +  6  rows.  Index  I  refers  to  the  production  industry.  Index  J 

refers  to  the  commodity  used  by  industry  I.  DATA  is  the  amount  of  com¬ 

modity  J  used  by  industry  I  in  the  production  of  its  product. 

The  indices  I  and  J  need  not  be”  integral.  For  example,  the  Bureau 
of  Economic  Analysis  classifies  the  United  States  economy  into  468  sectors, 
but  only  86  primary  industries.  For  this  data  base  there  are  468  values 
of  I  ranging  for  1.0  to  86.9999  and  474  values  of  J. 

Six  values  of  J  (0  to  -5)  are  added  to  t.his  page.  These  rows  are 

defined  as: 
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Figure  B-l.  Description  of  DSA  Data  Format 


J  Interpretation 

0  Undifferentiated  Imports 

-1  Wages  to  Labor 

-2  Burden  (Taxes) 

-3  Rental  on  Capital  (Property-Type  Income) 

-4  Capital  Inventory 

-5  Depreciation  (Total) 

The  term  "undifferentiated  imports"  refers  to  imports  which  have  no 
domestic  counterpart  or  which  are  consumed  by  final  users  in  the  form 
in  which  they  are  imported. 

The  sum  of  each  column  of  page  1  down  through  J  =  -3  is  the  total 
output  of  that  industry.  It  must  equal  the  total  use  for  that  product, 
i.e.,  it  must  be  equal  to  the  sum  of  the  corresponding  row  across  pages 
1-4.  The  data  for  J  =  -4  and  J  =  -5  are  not  required.  They  are  used 
only  if  the  user  wants  the  final  aggregated  data  put  into  model  format 
and  even  then  are  not  required  (see  Sec.  B.2). 


B.1.2  PAGE  2  -  INVESTMENT 

The  second  page  of  the  file  is  the  investment,  or  capital  formation 
page.  The  N  values  of  I  refer  to  the  industry  which  uses  the  capital. 

The  index  J  refers  to  the  commodity  required  in  the  manufacture  of  the 
capital.  No  labor  or  capital  is  required  in  the  manufacture  of  the 
capital  since  this  matrix  merely  describes  the  mix  of  items  comprising 
the  capital  for  each  industry.  One  extra  row,  J"=  0  for  undifferentiated 
imports,  is  provided. 


8.1.3  PAGE  3  -  CONSUMPTION 

Consumption  by  final  users  is  represented  by  page  3.  There  are 
six  final  users  to  which  the  data  may  be  applied.  The  user  may  define 
these  six  in  any  way  he  wishes,  but  the  aggregator  normally  assumes  the 
following  categories: 


* 
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I_  Final  User 

1  General  Population 

2  Military 

3  Federal  Government 

4  Local  Government 

5  Inventory  Change 

6  Undifferentiated 

If  the  user  wishes  to  add  more  consuming  sectors  he  will  have  to  modify 
the  model.  In  general,  however,  these  six  have  been  sufficient  for  all 
data  bases  used  to  date.  The  rows  of  page  3  are  once  again  the  industry 
codes  plus  J  =  0  for  undifferentiated  imports. 

B. 1.4  PAGE  4  -  TRADE 

The  trade  matrix  contains  just  two  columns.  I  =  1  refers  to 
exports  while  I  =  2  refers  to  imports.  The  DATA  value  for  J  =  0, 

I  =  2  must  in  this  case  equal  the  sum  of  the  undifferentiated  imports 
row,  i.e.,  all  of  the  undifferentiated  imports  are  "imported"  at  I  -  2, 

J  *  0,  K  =  4. 

B. 1.5  PAGE  5  -  DOMESTIC  TRANSFERS 

Some  data  bases  have  a  quantity  called  domestic  transfers.  When 
a  particular  industry  produces  a  completely  secondary  product,  the 
secondary  product  is  "sold"  to  the  primary  industry  for  that  product, 
which  distributes  it.  This  fictitious  transaction  is  called  a  domestic 
transfer.  After  the  data  is  aggregated,  the  diagonal  elements  of  this 
matrix  must  be  subtracted  from  the  diagonal  elements  of  the  production 
matrix,  page  1,  since  these  elements  represent  fictitious  transactions 
between  an  aggregated  industry  and  itself. 

B.2  USING  THE  AGGREGATOR 

The  aggregator  is  primarily  self-explanatory,  i.e.,  the  user  need 
only  answer  the  questions  asked  by  the  aggregator.  The  first  question 
asked  is: 
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HAS  DATA  BEEN  AGGREGATED  SUFFICIENTLY? 

The  aggregator  has  two  uses:  (1)  to  aggregate  economic  sectors,  and 
(2)  to  put  the  aggregated  data  into  model  format.  If  the  data  were 
previously  aggregated  and  put  onto  a  file  and  the  user  just  wants  to 
use  it  in  the  model,  or  to  change  it  slightly,  he  may  answer  yes.  If 
the  data  are  to  be  further  agcregated,  he  should  answer  no. 

The  aggregator  next  asks  the  user  for  the  name  of  the  file  which 
contains  the  data.  If  the  data  has  been  sufficiently  aggregated,  it  is 
read  in  and  the  program  informs  the  user  how  many  industries  are  con¬ 
tained  fn  the  data.  If  it  has  not,  the  program  needs  more  information. 
In  particular,  the  user  will  be  asked  to  supply  the  following: 

1.  The  number  of  industries  which  will  comprise  the 
aggregated  data  base. 

2.  The  industry  codes  which  go  into  each  aggregated  sector. 

No  industry  code  can  be  in  more  than  one  sector.  If  this  happens  the 
aggregator  will  inform  the  user.  The  codes  are  input  as  a  series  of 
code  ranges.  An  illustrative  sequence  is: 

Example  1: 

ECONOMIC  MODEL  DATA  PREPARATION— HAS  DATA  BEEN  AGGREGATED 
SUFFICIENTLY?  (Y  OR  N) 

yes 

INPUT  FILE  NAME 
f i 1 ename 

BEGINNING  TO  READ  IN  DATA: 

FINISHED  PAGE  1 
FINISHED  PAGE  2 
FINISHED  PAGE  3 
FINISHED  PAGE  4 
FINISHED  PAGE  5 

DATA  CONTAINS  (number)  PRODUCTION  INDUSTRIES 
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Example  2: 

ECONOMIC  MOOEL  DATA  PREPARATION--HAS  DATA  BEEN  AGGREGATED 
SUFFICIENTLY?  (Y  OR  N) 
no 

INPUT  FILE  NAME 
filename 

INPUT  NUMBER  OF  RESOURCES  (MAX  OF  100) 

10 

DO  YOU  WISH  DATA  TO  BE  INTEGRALLY  AGGREGATED  (Y  OR  N) 

HY  (ES)“  WILL  PUT  DATA  INTO  INTEGRAL  CODES  1-10  by  INTEGERIZING 
ALL  INDUSTRY  CODES1 
no 

INPUT  BEGINNING  AND  ENDING  CODES  FOR  INDUSTRY  1 
1  2.999 

ANY  OTHER  SEGMENTS  FOR  INDUSTRY  1?  (Y  OR  N) 
yes 

INPUT  BEGINNING  AND  ENDING  CODES  FOR  INDUSTRY  1 
3  3.9999 

ANY  OTHER  SEGMENTS  FOR  INDUSTRY  1?  (Y  OR  N) 
no 

INPUT  BEGINNING  AND  ENDING  CODES  FOR  INDUSTRY  2 
4.0  5.4001 

ANY  OTHER  SEGMENTS  FOR  INDUSTRY  2?  (Y  OR  N) 
no 


ANY  OTHER  SEGMENTS  FOR  INDUSTRY  10?  (Y  OR  N) 
no 


1 


This  is  useful  for  data  bases  containing  a  large  number  of  sectors 
divided  into  primary  industries,  such  as  the  BEA  US  data  base.  Instead 
of  the  user  having  to  input  code  ranges  1.0  -  1.9999,  2.0  -  2.999,  etc., 
the  model  automatically  sets  the  code  ranges.  This  will  aggregate  the 
data  into  the  orimary  industries  which  can  be  further  aggregated  later. 
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BEGINNING  AGGREGATION: 

FINISHED  PAGE  1 
FINISHED  PAGE  2 
FINISHED  PAGE  3 
FINISHED  PAGE  4 
FINISHED  PAGE  5 

After  the  aggregation,  the  user  has  the  option  of  listing  any  of 
the  aggregated  data  by  page  number,  and  of  changing  any  data  element. 

After  the  user  is  satisfied  with  the  data,  the  program  tests  for  consis¬ 
tency  between  row  and  column  sums.  If  inconsistencies  are  found,  the 
program  informs  the  user  of  the  problem  industries  and  allows  the  user 
to  look  at  and  change  the  data  before  testing  again.  If  the  data  is 
satisfactory,  the  user  may  copy  it  to  a  separate  file.  This  aggregated 
data  will  have  integral  industry  codes  as  opposed  to  the  primary  data 
base  which  need  not  have  integral  codes. 

B. 3  MODEL  FORMAT 

After  aggregation,  the  user  may  decide  to  have  the  data  put  into 
model  format.  Assuming  he  answers  the  question 

WOULD  YOU  LIKE  DATA  PUT  INTO  MODEL  FORMAT?  (Y  OR  N) 

with  "yes,"  the  aggregator  will  inquire  as  to  his  preferences  regarding 
the  number  of  capital  types,  the  names  of  consumables,  the  number  and 
definition  of  consumption  activities,  and  primary  data  arrays. 

B.3.1  CAPITAL  TYPES 

The  following  question  is  asked: 

CAPITAL  MAY  BE  DIVIDED  INTO  TWO  TYPES  BY  INDUSTRY  GROUP  IN  ORDER 
TO  DISTINGUISH  TWO  DEPRECIATION  RATES  (E.G.,  FOR  STRUCTURES  VS.  EQUIPMENT) 
—  DO  YOU  WISH  THIS  DIVISION?  (Y  OR  N). 


If  the  user  answers  with  "yes,"  he  will  be  asked  to  supply  the 
aggregated  industries  contained  in  each  type  of  capital.  Some  industries 
may  use  only  one  type  of  capital.  The  aggregator  will  determine  if  this 


is  so  and  delete  any  capital  not  required.  As  only  two  types  of  capital 
are  allowed,  the  total  number  of  capital  items  will  be  between  N  +  1  and 
2N  +  1  where  N  is  the  number  of  aggregated  industries.  Every  industry 
must  use  at  least  one  capital  type.  If  an  industry  is  found  which  does 
not  use  any  capital,  the  model  will  inform  the  user  and  abort. 

B.3.2  INDUSTRY  NAMES 

The  user  has  the  option  of  letting  the  aggregator  provide  names 
for  the  aggregated  industries  (IND  1,  IND  2,  etc.)  or  naming  them  him¬ 
self.  The  industry  name  may  not  be  greater  than  15  characters. 

B.3.3  CONSUMPTION  ACTIVITIES 

The  user  may  break  page  3  into  a  set  of  consumption  activities  in 
any  way  desired.  The  most  general  way  to  define  a  consumption  activity 
is  by  giving  the  industries  and  consumer  groups  (see  Sec.  B.1.3)  included 
in  the  activity: 

INPUT  GROUP  NUMBER  WHICH  IS  INCLUDED  IN  CONSUMPTION  ACTIVITY 

(number);  "S"  TO  STOP 

1 

INPUT  INDUSTRY  NUMBER  OR  RANGE  OF  INDUSTRY  NUMBERS  INCLUDED  IN 
GROUP  1 ' S  PARTICIPATION  IN  ACTIVITY  (number) 

5 

ANY  MORE  INDUSTRIES  FOR  GROUP  1  IN  CONSUMPTION  ACTIVITY  (number)? 

no 

ANY  MORE  GROUPS  IN  CONSUMPTION  ACTIVITY  (number)? 

no 

INPUT  NAME  FOR  THIS  ACTIVITY 

Contiguous  groups  may  be  listed  together'.  For  example,  if  the 
above  activity  were  to  be  total  consumption  of  industry  5,  the  user 
could  have  specified  groups  1-6  in  the  first  question  by  typing  in 
"1  6”  or  "1,6.”  It  is  not  necessary  to  list  each  group  separately. 

Conversely,  if  the  user  would  like  to  break  out  a  single  group  or 
combination  of  groups  for  each  of  a  list  of  industries,  he  may  answer  the 
second  question  as  follows: 
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INPUT  INDUSTRY  NUMBER  OR  RANGE  OF  INDUSTRY  NUMBERS  INCLUDED  IN 
GROUP  l'S  PARTICIPATION  IN  ACTIVITY  (number) 

5  10  I 

The  "I"  informs  the  model  that  industries  5  through  10  constitute  six 
independent  consumption  activities  by  the  same  group  or  groups. 

B.3.4  MODEL  PARAMETERS 

Besides  the  activity  levels  and  direct  requirements  matrix,  the 
economic  model  requires  the  following  data: 

(1)  Capital  inventory  levels 

(2)  Capital  depreciation  rates 

(3)  Current  level  of  inventory  growth 

(4)  Current  value  trends  (used  in  calculating  the  discount 

rates— not  used  directly) 

(5)  Discount  rates 

(6)  Gestation  times 

(7)  Population  growth 

The  aggregator  considers  which  of  the  above  it  can  calculate  and 
which  it  needs  to  ask  for.  For  example,  if  rows  -4  and  -5  of  page  1 
contain  data,  then  the  total  capital  inventory  and  overall  depreciation 
rate  is  known.  If  the  industry  contains  only  one  type  of  capital,  its 
inventory  and  depreciation  is  therefure  also  known.  If  two  types  of 
capital  are  used,  the  model  will  split  the  total  inventory  according  to 
how  the  user  wants  to  split  the  deoreciation  rates  of  the  two  types. 

The  aggregator  will  list  the  above  data  for  the  user,  inserting 
default  values  where  necessary  (e.g.,  gestation  times  are  defaulted  to 
one  year).  It  will  then  ask  the  user  to  change  any  of  the  data  he 
wishes  /ith  the  exception  of  capital  inventory  and  activity  levels. 

Some  of  the  data,  e.g.,  discount  rates  or  inventory  when  two  types  of 
capital  are  used  or  when  row  -4  is  blank,  are  listed  as  question  marks. 
These  data  will  be  calculated  by  the  aggregator  after  the  user  has 
specified  the  remaining  parameters.  Value  trends  are  defaulted  to  the 
growth  rate.  If  the  growth  rate  is  to  be  calculated,  i.e.,  it  is  listed 


I 

as  a  question  mark,  the  value  trend  will  be  listed  as  zero.  The  actual 
value  used  in  the  calculation  of  the  discounts  will  be  either  the  growth 
rate,  which  is  calculated  first,  or  the  value  input  by  the  user. 

After  the  user  has  changed  the  data  as  he  wishes,  the  aggregator 
will  calculate  the  remaining  data  (previously  listed  as  "?")  and  will 
relist  the  updated  data  arrays.  The  data  are  then  output  onto  file 
ECO. DAT  for  later  use  by  the  economic  model. 


APPENDIX  C 

GLOSSARY  OF  VARIABLES  USED  IN  DSA  ECONOMIC  RECOVERY  MODEL 
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ACON 


ACTLAB 

ARENT 


B 


CAPLAB 

CONLAB 

D 

DISCNT 

DLAMB 

DTAU 

EFFCAP 

EFFCST 

effLab 

EPSO 

FLAM 


Direct  requirements  matrix  —  ACON(I.J)  is  the 
amount  of  commodity  I  used  per  unit  of  activity 
J 

Labels  (names)  for  activities 

Amount  spent  on  capital  and  labor  per  unit  of 
activity  for  the  production  of  consumables  in 
the  current  economy  (first  NCAPM1  is  amount 
spent  on  capital  (each  capital  is  used  for 
only  one  activity  specified  in  ITRADE);  next 
NCON  is  amount  spent  on  labor) 

Exponent  for  utility  and  translog  functions 
(first  NCAPM1  for  translog  exponent  of  capital/ 
capital  trade-offs;  next  NCON  for  translog 
exponent  of  capital /labor  trade-offs,  next 
NCONAC  for  utility  exponent) 

Labels  (names)  for  capital  resources 

Labels  (Names)  for  consumables 

Depreciation  for  the  NCAPM1  capital  resources 
(D(NCAP)  *  population  growth) 

Discount  rate  for  the  NCAPM1  capital  resources 

Rental  value  for  capital  resources  (each  time 
period) 

Difference  between  TPER  and  TLAG  (DTAU(J)  = 
TPER*(1 .0  -  AM0D(TLAG(J)/TPER,1.0) )  ) 

Effective  capital  inventory  (used  when  two  types 
of  capital  are  Used  to  make  one  consumable) 

Effective  cost  of  producing  capital 

Labor  used  in  production  of  particular  consumable 
(not  constant  but  varies  with  time  period) 

Epsilon  parameters  used  .in  Everett  algorithm  for 
phase  two 

Shadow  value  (offset  by  Dne  time  period) 
(FLAM(K=1))  contains  horizon  DLAMB  values  for 
capital 
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FUDGE 

GAMMA 


GROWTH 

ITRADE 

K 

NACT 

NCAP 

NCAPM1 

NCON 

NCONAC 

NCST 

NPER 

NPRDAC 

PERLAB 

PRMLAM 

R 

REQ 


Normalization  coefficient  for  translog  function 
(first  NCAPM1  for  capital /capital  trade-offs; 
next  NCON  for  capital/labor  trade-offs) 

Coefficient  for  utility  and  translog  functions 
(first  NCAPM1  for  translog  function  of  capital/ 
capital  trade-offs,  next  NCON  for  translog  function 
of  capital/labor  trade-offs  (if  GAMMA  .GT.  0  IT 
IS  GK  —  if  .LT.  0  ABS( GAMMA)  is  GL,  GK  +  GL  =  1.0); 
next  NCONAC  are  utility  coefficients  (activity  = 
GAMMA/ (cost**B0))  ) 

Fractional  growth  of  industries  producing  capital 
in  current  economy 

Array  stating  for  which  industry  capital  J  is 
used 

General  index  for  time  period 
Number  of  activities 

Number  of  capital  resources  (including  labor) 

(NCAP  must  be  .GT.  NCON  .AND.  .LE.  <>*NC0N  +  1) 

Number  of  capital  resources  (excluding  labor) 

=  NCAP  -  1 

Number  of  consumable  resources 
Number  of  consumption  activities 
Total  number  of  resources  *  NCAP  +  NCON 
Number  of  time  periods 

Number  of  production  activities  =  NCON  +  NCAP  -1 

Labels  used  when  printing  out  arrays  with  time 
period  indices 

Lambda  prime  -  gross  profit  in  the  manufacture 
of  commodity  I  which  is  distributed  between 
capital  e  labor 

Resource  level  for  capital 

Equilibrium  levels  for  capital 
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RFUT 

RLAST 

RORGNL 

RSCALE 

RTRGT 

SCRTCH 

SUMLAB 

SUMMRY 

T 

TCONS 

TLAG 

TPER 

W 

X 

z 


Resource  level  one  gestation  time  in  the  future 
if  no  investment  activity  occurs  (used  in  phase 
one  only) 

Resource  level  achieved  in  previous  pass  for 
given  time  period  and  capital  item 

Starting  level  of  resources 

Scaling  factor  used  in  phase  one  method  of 
determining  investment  activity  levels 

Target  resource  level  for  one  gestation  time 
in  the  future  (used  in  phase  one  only) 

Scratch  pad  array  —  used  for  a  variety  of 
functions 

Labels  used  in  summary  tables  at  end  of  run 

Summary  tables  array 

Time  at  center  of  time  period  K 

Total  consumption  of  commodity  I 

Gestation  time  for  capital  (years) 

Length  of  time  period  in  years 

Storage  array  which  contains  every  other  array’, 

W  is  written  onto  insurance  file  after  each  pass 
for  restarts 

Activity  level 

Minimum  consumption  activity  level 


94 


APPENDIX  D 
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